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Abstract 

We find general relations between RG equations and planar unitarity-analyticity. These relations are summarized in 
meromorphization procedure, generalizing the Pade approximation in the limit of infinite order. 

We also investigate confinement conditions for the mass spectrum in asymptoticaly free QFT and lay down systematic 
framework for a expansion suggested in previous papers. The new relations for meromorphization of symmetric conformal 
tensors are found, with resulting rich mass spectrum as a function of spin n. 

Explicit intergal representation for the triple string vertex F is found. This corresponds to resonanse theory with infinite 
number of masses and Lagrangean 
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1 Introduction. 

Recently U , we revived old approach to planar confining QFT [31 0] . The motivation for this revival was the observation [8 that 
the Ads^ / C FT theory being regularized by a sharp cutoff in the vicinity of its -R4 boundary leads to the mass spectrum exactly the 
same as suggested in [SJIJ, namely roots of Bessel functions. It was proven only for conserved currents 2-point functions of CFT, 
where there were no corrections to anomalous dimensions of conformal operators. However, the regularization in [31 [3] following 
solely from unitarity and analyticity, there are strong reasons to believe that general formulas, with anomalous dimensions as 
indexes of Bessel functions also apply to regularized Ads /CFT theory. This would be an exciting phenomenon: exact relation 
which holds in (regularized) Ads /CFT theory to all orders in coupling constant. Recent nonperturbative computations [9] of 
anomalous dimensions in A'^ = 4 SYM theory as functions of coupling constant make this conjectured relation even more exciting. 

Can we now do something about the mass spectrum in confining QFTl At phenomenological level we know that Pade 
method works for QCD with vacuum condensates (016]). 

Before presenting the results of this paper, let me explain how and why the mass spectrum is related to the OPE. It is 
relatively simple to go from the spectrum to the OPE, but not so simple to go the other direction. 

Take the sum of pole terms 




(1) 



and apply the Euler summation formula 



1 - exp((}„) [m2 - t 
This produces asymptotic expansion (with Bp being Bernoulli numbers) 



(2) 



J n=0 




(3) 
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with the first term corresponding to replacement of discrete sum by an integral, and the rest being corrections. We can 
introduce the density of the spectrum, assuming we know analytic continuation of Z„ and m^j to continuous values of n 

p{m') = -^,Z{ml)^Zn, (4) 
and rewrite this expansion as an integral plus corrections 

The first term would produce the continuum spectrum, corresponding to leading order of asymptotic freedom. The higher 
Bernoulli terms will produce powerlike corrections, proportional to negative integer powers of t. Those must correspond to the 
higher terms of OPE, or else they must vanish. Take the simplest possible Anzatz, the notorious psi function with 

Z^ = l,ml = n + l (6) 

It corresponds to constant spectral density, and so, the only corrections to the leading term Ini will come out as integer 
negative powers. That clearly contradicts the spectrum of OPE for large N YM, where there are no conserved currents except 
Energy Tensor. All dimensions are known to get perturbatively renormalized, so that there must be calculable fractional powers 
of Ini in front of negative integer powers. 

Such terms can only come from the integral term, due to some nontrivial behavior of discontinuity 

5G{t) ^ Trp{t)Z{t) (7) 

which must posses all these powers of t and Int. At the same time we observe that all the Bernoulli terms must vanish 
identically, otherwise there would be fake operators with integer dimensions in our theory. This provides highly nontrivial 
restrictions on the mass spectrum m„ and residues Z„. 

Thus, we see that there are some hidden connections between the spectrum of anomalous dimensions of perturbation theory 
and the physical mass spectrum, and given the mass spectrum it is straightforward to compute the UV asymptotic expansion 
and compare it with the OPE, and thus verify correctness of the hypothesis about spectrum. 

Our goal in this paper is to elaborate the same hidden connection between masses and anomalous dimensions going the other 
way- from asymptotic freedom to the physical spectrum. This way is much harder, but not impossible, as I have shown 35 years 
ago. 

In this paper we continue investigation of Pade regularization. We rederive, debug and reinterpret old formulas of l3i|4j, and 
we make some new advances and insights. We are not going to assume conformal symmetry but rather consider an arbitrary 



Z{s) 



(5) 
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large N QFT with confinement, i.e. discrete spectrum of states in every channel. The specific structure of CFT and its Ads 
correspondence will not be used here. 

As it was discussed at length in [3l 1^ one must use the matrix Fade approximant rather than scalar one. Nothing changes 
in principle except necessity to keep track of order of matrix multiplication. We carry out this matrix Fade approximation 
to the end for the case of arbitrary symmetric conformal tensors, where the 2-point function is the matrix in space of 0{d) 
representation. As a result we get new rich mass spectrum for such conformal tensors, overlooked in previous work. 

The essence of our method is to impose on the n-point functions in momentum space their correct analytic properties, i.e. 
meromorphicity. In case of 2-point function G{t), where t = —p^ we demand that it only has poles, all located on the right 
semi-axis. We shall suppress the matrix indexes so far and restore them later. We call this transformation meromorphization: 

M[G,Q,t] = ^, (8) 

P{t) = G{t)Q{t) - / 5,G{s)Q{s). (9) 

where dsG{s) is discontinuity of G{t) across the cut at positive real axis s > 0. The entire function Q{t) here is chosen in 
such a way that the dispersion integral here decreases faster than any power, i.e. all the Laurent expansion coefficients of this 
integral at infinity must vanish: 

/>oo 

ds5sG(s)Q(s)s" = 0,n = 0,1,... (10) 





The motivation for this requirement is that these negative power terms in asymptotic expansion of G{t) come from the OPE, 
with coefficients proportional to VEV of various operators with vacuum quantum numbers, such as traces of powers of F^j, in 
YM theory. We would like to preserve all such terms while adding our corrections, so that these corrections are exponentially 
small at large momenta. Note, that in case (which we would like to eventually achieve) when the mass spectrum is given by the 
roots of Q{t) the extra term we add will become identically zero, because 5sG{s)Q{s) cx S{s — m^)Q{s) = in this case. 

In other words, meromorphization is an identical transformation of the full theory, with the purpose to improve 
the perturbation expansion. The physical meaning of meromorphization is to declare that the theory is an infinite collection of 
free fields $j;(x) with various masses and couplings to the gauge invariant YM sources Jk{x) conjugate to operators Ok{x) in 
YM theory. 

5e//= J d^xJ2'^aMabiiV)<i>b + 'i>afab{^)Jb (H) 

ab 

The mass spectrum is given by vanishing eigenvalues of infinite matrix M(p) and the 2-point function will factorize as 
r^{p)M~^{p)r{p). As it was discussed in the previous paper, the Matrix Fade approximation of G provides precisely the same 
representation, disguised in a form of right or left multiplication 

G = f'iM-^f^PQ-^ ^Q-^P. (12) 
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The asymptotic results at large order M, N of the approximant are the same, and they can be proven to be factorized as required 
by free field theory, with real masses and real coupling constants. This remarkable coincidence follows from the general theorems 
of Pade theory for the Stiltjes matrix functions. In particular, positivity and factorization was explicitly demonstrated in the 
previous paper. There are more subtle theorems, such as monotonous decrease of approximated mass spectrum with the rank 
N of Pade approximant, and theorems about absolute convergence of approximant in complex plane of t. 

Do not let the word "approximation" mislead you. The so called approximant is a rational matrix function, with coefficients, 
satisfying certain linear matrix equations: so called Pade equations. These equations are usually solved numerically, which we 
are not going to do. We rather find exact solutions of these equations for arbitrary order M, N of numerator and denominator 
polynomials. This solution represents certain integral transformation of the matrix function G{t) which transformation can be 
defined exactly in every order of perturbation theory of OPE. Then we find drastic simplification in the limit of large M, N 
which allow us to go further with the mass spectrum computation. 

We start from Pade approximation of G{t) near large Euclidean point t = —A in the limit when the orders M, N of Pade 
approximant both go to oo at fixed = 2AIN ^ ^j^jg jj^s dimension of length, so that conformal symmetry is explicitly 
broken. It may be restored in the limit i? — oo . In case of precisely conformal field theory, such as SYM dual to Ads string, 
we are looking for dynamic mechanisms enhancing this conformal symmetry breaking (by VEV of composite fields) so that it 
persists in a limit R ^ oo . This is similar to introduction of small magnetic field into ferromagnetic below Curie point to obtain 
spontaneous magnetization in a limit of vanishing magnetic field. 

In case of running coupling constant under consideration in this paper, the conformal symmetry is broken already in the 
second loop. It would be convenient, however, to separate the effects of renormalization of anomalous dimensions, preserving 
conformal symmetry of YM theory, and effects from running coupling constant leading to asymptotic freedom and - as we all 
hope - confinement. We study all these effects in some detail in this paper, to review and extend the old results ^(4^. 

Few words about relation of two large space scales: confinement scale Rqcd = which is exponentially large in 

perturbation theory, and our regulator scale R which we eventually must tend to infinity to recover original theory. We start 
with the opposite limit of i? ^ Rqcd which is like placing our QFT in a box size R much smaller than confinement radius. In 
this limit the spectrum is obviously discrete in every order of perturbation expansion in running coupling A^^ ~ 1/ ln{RQCD/R) 
. Then, we start increasing R and develop the technique of extrapolating the mass spectrum to the physical limit of i? ^ Rqcd 
in which limit we recover original theory. 

The general theorems of Pade theory for Stil'tjes functions are crucial for this extrapolation. These theorems state that the 
position of poles of approximant monotonously decrease as functions of its order. In our case this means that mass spectrum is 
monotonously decreasing as R passes the confinement scale Rqcd and goes to infinity. Nothing pathological like Landau pole 
can happen with this definition of running scale. Confinement would correspond to a finite limit of every mass in the spectrum 
at R = oo. We find the way to redefine the running coupling a so that it also reaches the finite limit, namely a = 1. We argue 
that series expansion in this running coupling a has finite radius of convergence, (precisely a = 1). The expansion terms are 
calculable, so that the problem is to extrapolate the series to the weak singularity (infinite first derivative) at the convergence 
radius. 

This approach was initiated in the old work ([31 H]) but here we find new powerful method to compute terms of expansion. 
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and also correct some old errors. 

We derive all the relevant meromorphization formulas from scratch in Appendix, correcting some errors and typos along the 



way (including some which propagated into the recent paper W ). The end results in the limit M, N oo, 
are quite simple. 



2MN 



~ const 



First, let us study the relation between P{t) and Q{t). The two terms can be combined as a single contour integral 

dsG{s)Q{s) 



Ct 



2m{s - t) 



(13) 



where Ct encircles anti-clockwise both the pole at s = t and singularities of G{s) at positive real axis. 

The contour cannot be closed in the left semiplane because of exponential growth of Q{s). The collision of singularities at 
t — ^ in this contour integral does not lead to any singularities because they do not pinch the contour. So, P{t) is an entire 
function as well. Its expansion coefficients p„ in powers of {—t) are given by the following contour integral (in all subsequent 
formulas we choose i? as a unit of length) 

dsG{s)Q{s) 



Pti 



Co 



2TTi 



-(-.)- 



(14) 



6 



Let assume that we know the Melhn transform of G{s) 

G{t)^j ^^F[u){-tr. (15) 



where Cp encircles the singularities of F{uj). Substituting this into relation for p„ we get 



The last integral here is nothing but an expansion coefficient qn of Q{t) analytically continued to complex values of n 

oo oo 

= £ ^^H9m— (19) 

This is quite a remarkable relation. Once the coefficients g„ of the entire function Q{t) are known as analytic functions of 
index n the other entire function P{t) is as good as known. We shall see below how this relation works in practice, within planar 
graph expansion. 

In general case, of many point function, it can be represented within planar graph expansion as multiple Mellin integral |10] 
of various planar kinematical invariants 

Sab = - ' (20) 

(s.) = / Dn{Lo)F'' {uj.) W [SabT-' . (21) 

<ab> 

Meromorphization in variable Sab would correspond to multiplication by Q{sab) and contour integration as before. We will 
get multiple power series expansion in numerator and product of Q in denominator 

Dn{u)F^ (CO..) n , (22) 

Below, we carry out this procedure to the end for conformal 3-point vertex. 
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2 Conformal Approximation. 



The equation for g„ is less trivial than equation for p„. Let us assume that F{co) is given by sum of two terms 

F{Lo) = ^—+rS{uj), (23) 
w — 7 

where 7 is the extreme right pole of F{lo) with the residue r. In terms of G{t) this corresponds to powerlike behavior. The 
second term (which can also have singularities at w = 7) will be treated as perturbation. This term, to be studied in the next 
Chapter, sums up all the effects of running coupling constant, and starts from the second order in perturbation expansion. As 
it follows from general relation derived in Appendix , the coefficients satisfy the following equation 

^ 00 

'^^ " r(M + i)r(Ai + 7 + i) ~ ^ ' ^^^^ 



n=0 



f duj 

Pfi = rq^—i 2^i^^^'^'^^~" ^^^^ 

Sij.n=(f ^F'^{n,n + u-^)S{ijj), (26) 
Jcs 

sin(7r6) sin(7r (6 + 7)) T{b + l)T{b + 7) 
7rsin(7r7)(6- a) T{a + l)T{a + ^) 

The same formulas also provide analytic continuation for complex values of fj, . 

Let us consider first the leading term, corresponding to conformal theory, and let us now restore the length scale R 



^° r(/i + i)r(M + 7 + i)'^° r(M-7 + i)r(/i + i)' ^^^^ 

gO(f) = {-t)-il-y (2V-tR^^ ,P°{t) = r{-t)il-j (2V -tR'^) ■ (29) 
The meromorphized 2-point function reduces to the ratio of Bessel functions 



^ ' i-tyilj {2^f=tW) V 7^/7 {2V^tR^) J ^ ' 
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In special case of integer 7 — )■ n corresponding to conserved currents, the residue r must grow inversely proportional to 
sin(7r7) 

r' 



r = 



sin(7r7) 

so that (up to irrelevant additive regular term r(— t)" ) the 2-point function becomes 



U ' ' t ' /„(2V^l7P) 



with first term corresponding to its conformal limit. 



The physical limit corresponds to large R where the correction to power term decays exponentially, as exp (— 4\/— ti?^), 
except at positive t where G^{t) has infinite number of poles with positive residues. These properties follow from the general 
theorems of Pade theory for Stiltjes functions (for readers convenience the proof of positivity was reproduced in ). The 
physical meaning is transparent: these are the composite states of our large N theory in conformal approximation. 

Note that exponentially decaying corrections to the scaling limit in momentum space correspond to analytic terms in coor- 
dinate space, going in powers of ^. Such corrections do not correspond to any physical operators in OPE of CFT. From the 
point of view of confining QFT with conformal symmetry explicitly broken by the beta function we do not see any problem with 
analytic terms like these. In the limit of i? — 00 these terms decay as negative powers of R , but after the summation of a 
perturbation expansion (see below) these terms may stay finite. 

3 Perturbation Expansion for Mass Spectrum in Dimensional Regularization. 

It is most convenient to study the running coupling constant with dimensional regularization scheme. In that regularization, 
the 't Hooft coupling constant A — Ng^ has dimension of mass^^ where (in our notations) e — 2 — d/2. The 2-point function 
expands in power series in A with coefficients being some functions of e times powers of momentum as dictated by dimensional 
counting. Ignoring spins, we have 

00 

G(t) = ^A^-Vfc(e)t^'^ (31) 

The term ft represents the sum of fc— loop planar diagrams for the 2-point function. This corresponds to Mellin transform 
F(uj) being a simple sum of pole terms 

F(c.) = y (32) 

k—l 
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which correspond to r = /i(e), 7 = 



— e in above equations. The equation for denominator becomes 



^ 00 

= r(M + i)r(M-e + i) - § ^^^^ 

V = E %^r-.(M,n-me)(Ai?2^)'". (34) 

m=l •'1'^^'' 

We can write down recurrent equations for expansion coefficients of in powers of the 't Hooft bare couphng. 

00 



9/. 



(r) 



J2hl^J2r-e{f^,n-me)q^:-^ (37) 



m=l ■^l(^) n=0 



The expansion for the masses around the roots rrii of Bessel function can be done by iterating equation 



^^2^ng„ (i? mj = " 2^ ^ J (38) 

n=l n=0 

Note that (as it should have happened) the overall normalization of G{t) corresponding to multiplicative renormalization 
constants of the operators, dropped in our equation for the denominator. In particular, the one loop term /i(e) had the pole 
at e = related to the logarithmic divergence of the one loop integral in d = 4 dimensions. Now this pole enters denominator, 
so that, in effect, we have an extra factor of e in the numerator. This cancels the pole at e = 0, coming from the kernel 
F-£ (m, n — me). Namely, there is the factor sin(7r7) = — sin(7re) in denominator. Now this pole at e = is compensated by 
/i(e). So, the counting of e— poles remains the same as in momentum space, though, of course extra zeroth and positive powers 
of e coming from the kernel, change the; r(;sulting finite terms. 

As we learned in the seventies from the famous 't Hooft's work, in order for the observables to remain finite in the limit of 
e = all we need to do is to renormalize the bare coupling 

A = R-^'\r (1 + ci{e)\Ti + C2{e)\\ + ...) (39) 

The physical coupling A/j corresponds to space scale R and is supposed to remain finite at e = 0. For that, the coefficients 
Cfe(e) must have some poles, which are designed to cancel the multiple poles in Laurent expansion of coefficient functions 
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/fc(e) ^ e^'^ + ... + + 1, after some other multiplicative renormalization of G{t). The existence of such universal functions 
Cfe(e) that cancel poles in all 2-point functions is called renormalizability. We take this for granted here. The UV regularization, 
dimensional or otherwise should not be affected by the IR regularization such as ours. 

Moreover, we assume these functions known. Every multi-loop calculation of the YM theory by necessity produces these 
functions, and also the similar non-universal functions in the renormalization constants for the gauge invariant operators. 
Naturally, we choose our parameter i? as a mass scale. This will produce, after expansion of all the powers of t and cancellation 
of the poles in e in every order in Xn some linear combination of powers of In R'^t in the momentum function G{t). 

As for the expansion of our denominator Q{t), it will produce some universal numbers for expansion coefScients q"^ when 
re-expanded in running coupling A/j. The powers of R will all cancel, as we have chosen R as our physical scale in the definition 
of the running coupling constant. This follows now from trivial dimensional counting: there are no dimensional parameters left. 

Let us see how this expansion starts in the lowest order. We take the first approximation, by multiplying l^-e (a*, n — e) by 
q'^'' . We get the following calculable sum 



At integer fj, = m 



Comparing this with 



(1) ^ /2(e)sin(27re) ^ T{n-2e) 1 

7r/i(e)r(M + l)r(/x - e) r(n + 1) (n - e - m) ^ ^ 

^ /2(e) sin 7r(e - /i) 

/i(e)r(/x + l)r(/x + e+l)sin7r(e + M) ^ ' 

we see that the first term in qm^ leads to overall renormalization, but the second term, being combined with the same term 
in q^ is equivalent to effect of anomalous dimension 

Z = l + fM.„, ,45) 
7--e+-^A«^-^AH. (46) 
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Therefore, the masses to the first order are given by the roots of J^{mR). 

This is, of course, to be expected: the second loop adds one more logarithm with yf^^ in front. This is equivalent to 
anomalous dimension to this order. The running coupling constant displays itself in the next loop. The nontrivial fact is that 
we derived this from the Pade equations, not from the conventional techniques of summing up logarithms. 

Higher terms are straightforward to generate, though I do not know whether higher order sums would be analytically 
calculable, like in the first order. These higher order terms involve some hypergeometric sums, which can be expanded in powers 
of e producing higher PolyGamma functions Psi(n + 1, fc) . The good news, however, is the convergence factor we see in this 
perturbation expansion. At large m our factors decrease, at least they do so at any finite e. 

1 

F-e (/i,n-me) — 



The problem of growth of planar graphs of high order (so called renormalons) was not specific to four dimensions: at any 
finite e 't Hooft analytic argument would work as well. This argument does not apply to Pade regularization, precisely because 
it restores correct analytic properties of 2-point function. The 't Hooft 's condensing singularities in complex A plane correspond 
to these same poles we have here, but old perturbation expansion did not have. Nothing negative can be said about positions 
of these poles from the point of view of analyticity: the renormalon argument loses its ground. 

Another important comment. With the dimensional regularization scheme we adopted here, there is no need to separate 
effects of multiplicative and additive renormalizations. The 2-point function of composite fields has both of these effects, as 
it is clear already in the leading order. There is a logarithmic term, and logarithm transforms additively under rescaling of 
the momentum cutoff. Therefore, the CS equations should be used with some care. These equations do not apply to the 
2-point function in momentum space, but rather to coordinate space at non-coinsiding points. The ill-determined subtraction 
polynomials in momentum, which reflect the additive renormalization, drop after Fourier transformation. Another, more honest 
alternative would be to write the CZ equations for the discontinuity of the 2-point function across the cut at positive t. This 
also eliminates all the additive renormalizations. Discontinuity of the logarithm, for example, is equal to tt and does not change 
with rescaling of the cutoff. 

But in dimensional regularization there are no additive renormalizations to begin with. Every Loop produces just a power 
of momentum with well defined coefficient in front. This leads to great simplifications in meromorphization. As we have seen 
above, these power terms are directly translated into calculable terms in our Pade equations. With powers of logarithms, one 
would have to represent those as derivatives of powers to compute, which would lead us back to dimensional regularization. 



4 Confinement and Beta Function. 



Let us now discuss the confinement condition, which arises in the limit i? — > oo. According to the previous papers [HII] taking 
this limit of mass spectrum is equivalent to its minimization over R in view of monotonic decrease of masses as functions of 



12 



Pade order. One can write each mass as follows (in convenient logarithmic scale, in units of QCD mass squared scale I^qcd)'- 



In = lim (- In {R'Kqcd) + Pr{\R)) (47) 

where is the running coupling, satisfying the RG flow equation (negative sign because of space scale R instead of mass 
scale) 

and Pi{\R) is given by our regularized perturbation expansion, starting with constant: 

P.(Afl)-p.(0) + Aflp^(0) + ...; (49) 
P.(0)-lnr2 (50) 
Jo in) - 0; (51) 

^ -m^^' 

In order to compute these terms of expansion, one has to perturb the Bessel solution for Q by higher terms as explained in 
the previous Section. Also, one should expand the anomalous dimensions, entering the indexes of the Bessel functions, in power 
series in . Combining all these perturbations, we get the higher terms of expansion (|49l) . 

Replacing the limit i? — )• oo by the extremum condition (valid for our monotonous function!) we get 

In (^^\ = min (- HR^Kqcd) + P^{\R)) (53) 

K^QCDj R 

In order to make this limit non-singular we introduce Lagrange multiplier a as foUowfS 



^The transformation <3>(a;) = Ui.\ny{F{xj)—xJ{y)) is called Legendre transformation. It is well defined provided the extremum equation F'(y) = xj'{y) 
has only one solution, which is our case. Once this is true, there are several nice properties of Legendre transform. In particular, as the value of the 
function at its extremum does not depend upon the choice of the variable, the function <3?(x) is invariant with respect to the change of initial variable 
y : y G{y). In particular, this means invariance with respect to renormalization scheme. Another useful feature is the equation for derivative of 
: ^'{x) = —J{y*), where y*{x) is the position of extremum. The other terms in derivative y' {x){F' (y) — xj'{y)) are identicaly zero due to the 
extremum condition. In our case 

x = a^,.J{y) = \n{R^AQCD), 

and we can either use ln{R'^ Aqqj}) as independent variable y, in which case J{y) = y and F{y) = p{Afl), or else we can use y = Xji, in which case 
F{y) = p{y) and 

, , /■^« d\ 

Results are independent of the choice of variable. 
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In 



A, 



QCD 



Jnn_^ jniin {-a^ \tl{R^Kqcd) + Pi{\b)) } 



(54) 



Before we start working with this modified equation let us discuss in some detail the crucial point of monotonous decrease 
of masses with respect to R. For a = 1 it follows from the Fade theorems. However, at a < 1 this is no longer true. In case 
of confinement the function pi(Afl) must grow at large R precisely as ln(i?^) to cancel the first term. In other words, in case of 
confinement we have not only the inequality 

^/^^(^fi) < 1 f55) 
dln(i?2) ^ 1 ^^^> 

which reflects monotonous decrease with R of original mass m| in (1531) , but we rather have an equality 

dpi{\R) 



rfln(i?2) 

This means, by continuity, that at least at large enough R 



dHR^) 



(56) 



(57) 



so that the minimum is unique. We cannot prove this positivity condition for all R but we can prove that it holds in the opposite 
asymptotically free region of large negative \\\{R^Kqcd)- In this region the masses arc roots of Bessel functions J^{mR) — 0, 
with index 7 being anomalous dimension of the operator (up to A|j corrections). So, we have 



dpi{\R) 
filn(i?2) 



-2/3(A. 



-^/3"(0)P'(0) 



4/3"(0)7.:(0) ^ 



> 



(58) 



7=0 



The derivatives of roots of Bessel function with respect to its index at zero index are known to be positive, and so is — /3"(0). 
The anomalous dimensions for some non-conserved currents are positive in the first order, so that 7'(0) > . We can choose 
which root we take as a definition of our transformation from R to a. The parameter R is universal, so we can express it in 
terms of a from the lowest scalar mass (7'(0) > in that case). Then we take ratios of mi/rriQ and re-expand the Xr expansion 
in series in a. 

As a result, we see that is positive in asymptotically free region as well as in the confinement region. One could 

imagine some pathological behavior with this derivative changing sign two times along the way from asymptotic freedom to the 
confinement region. In that case, there may be several branches of the solution for Xr as a function of a. The best we can do 
is to compute the perturbative branch, available to us, as a series in a and check whether there are some phase transitions on 
the way from a = to a = 1 by estimating the radius of convergence. 



14 



Let us now dwell on the minimality condition, which can be rewritten in terms of effective coupling Xj^ 
and the equation for the logarithm of mass which can be rewritten as follows 



(59) 



In 



A 



QCD 



' dX 



(60) 



In order to eliminate the unphysical variable A r we must expand it in power series in a by inverting (|59|) and substitute into 
(pO)) . With the now popular choice of beta function, suggested in the seventies ( 4j, formula (1)) 



Aqcd = a exp - 



5 In Aq 



/3(A) 



a h 

"a2 + 3^' 

967r2 _ 102 

11 ' ~ m' 



we have 



In 



A, 



QCD 



(l 



61nA 



/'o(Ai?) 



(61) 
(62) 
(63) 

(64) 



Q = l 



It is important to note that after taking minimum with respect to R, already at arbitrary a < 1 the dimensional transmutation 
took place. The IR cutoff R disappeared, so that the mass has the correct QCD scale times some function of dimensionless 
parameter a 



4 = i^QCoFia) 



(65) 



These equations provide the basis for systematic expansion oi F{a) using the perturbative expansions for /3(A), 7(A) and po{X). 
The effective coupling starts linearly with a 



Xr = a 



+ 0(a2). 



F{a) = (^1 + 4a^ap',{0)^ + 0{a') 



(66) 
(67) 
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After that, the dimensionless rations mi/niQ can be expanded in a, starting with ratios of roots of Bessel functions Jq 



Inf— ) =^(A(AH)-po(Ai^)) (68) 

= In (^^) + - Pom + 0{a% (69) 

Jo(r.)=0; (70) 
d In r, 



pK0)=7K0) 



^7 



(71) 



7=0 



Since a starts hnearly in A we may view this a as a physical running couphng constant. The advantage of this redefinition of 
couphng constant is obvious: unhke original coupling A this a tends to 1 in the strong coupling limit R = oo. In the weak 
coupling limit a goes to zero as — 1/lni? , same as Xr. 

The reader may wonder: how did we get around the notorious Landau pole? The effective coupling A_r can have any 
singularities as a function of R, it can have a pole, or even a branchpoint such that one would not be able to continue Xr beyond 
some value of R without getting imaginary parl|3.These troubles do not reflect the physics of our system but rather the poor 
choice of effective coupling. We know from Pade theorems that In(mass) decrease monotonously as functions of Ini? so we can 
perform Legendre transformation and rely on nice properties of masses as functions of Legendre parameter a^. 

In fact we can plot this dependence qualitatively, as it was done in my old papers. This follows from the famous Legendre 
formula, in our case 

-^^-2aHR^A,cn)^-Aal^ M ^''^ 

At small a the integral is dominated by the lower end, and produces ^ which cancels a and leads to a positive finite 

limit. 

d In rriQ Aaa 
da Xr 

It then reaches the maximum at Aqcd 



AJap'oiO) >0;a^0. (73) 



— -^ = 0,i?2^- ; 74 

da Aqcd 

This is all what happens with our theory at would-be Landau pole. This is a position of the maximum of each mass in the 
spectrum as function of our physical coupling a. After this the mass goes down and reaches finite limit at a = 1. However, this 



This is what happens with the popular minimal choice of beta function at Xcrit - 
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0.2 



1.0 



Figure 1: Inm^ = a — {1 — a) ln(l — a) ; < a < 1 



is a singularity, as the derivative 

din m?, 

— ; = — oo; a ■■ 

da 



Here is the simplest function with such a behavior: 

The extrapolation to a singular value a = 1 can be done by the continuos fraction made from available few expansion terms 
(hopefully, 10 at modern level of analytic computations of the Feynman graphs of massless theory). The convergence of continuos 
fraction at a = 1 to exact value m = 1 for this prototype function goes as follows (for various orders of continuous fraction) 
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order 


In m(a = 1) 


2 


2 


3 


1.6 


4 


1.25 


5 


1.17391 


6 


1.11111 


7 


1.08485 


8 


1.0625 


9 


1.05049 


00 


1.0 



Let us also discuss the case of perturbations around the nontrivial fixed point A*. In this case the beta function vanishes 
Unearly so one has to use another parameter ^ = ^/a. 

^mdxto + ^to = 0, /3(A) = (A - A*)/?' + O ((A - A*)^) (75) 

and expand in series of ^. The rest of the arguments goes the same way, with replacement of — >■ ^. So, the CFT perturbed 
around the fixed point by a running coupling constant can confine as well, and the above expansion can be used to compute its 
spectrum in terms of expansion in powers of ^. 



5 Meromorphization of Conformal Vertex 

From the point of view of string theory (or, better to say, dual resonance theory, as we do not know nor we need to know explicit 
string model), the denominator Q provides quadratic part of the effective Lagrangian 

$]$aQa6(V^)$6. (76) 

ab 

We ignore here the tensor structure, which in momentum space depend of = In order to achieve analyticity, 

corresponding power of must be present in Q to cancel kinematical singularities of \p\ = ^f—l in denominator of . We 

shall not go into these details in present paper for clarity of presentation. We shall take here normalization of 2-point function 
as pure power of x in coordinate space, so that in momentum space (for scalar case) 

G{-e) = j d'^xx-^^ exp(ifcar) = a(t/)fc2-; a(!.) = f^dTV) ' = ^ ~ f " ^^^^ 
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Let us now introduce some cubic interaction 



J2 Tahc {-idi, -id2, -ids) $a(l)$b(2)$c(3) (78) 

abc 

and compare the conformal 3-point function 

{Oi{xi)Oj{x2)Ok{x3)) = ajfca^if *-^^+^'=X2-3^^-^''+^^X3-i^'=-^^+^^ (79) 

to the tree diagram of this theory. We need to go to momentum space, meromorphize and multiply by throe propagators Q 
to obtain the vertex Tatciqi, li, Q3) as a function of external momenta qi. Ignoring tensor structures and skipping indexes a, b, c 
we have triangle diagram with power propagators (here and below we denote by i,j,k the cyclic ordered indexes 1, 2, 3) 



G(si,S2,S3) = Ci23<7(^)a(^)a(^) J - q^' \q - q^\^^ \q - q^' -,3^ = ql,a, = A, + - \ - d. (80) 



Exponentiating these power propagators we get 



TM 

The Gaussian integral over k is straightforward. After some algebra (thanks to Mathematica ) we get simple expression 



This integral further simplifies by the change of variables from Xi to 



Diyi,y2,y3) „ „ . s/yrnvd. , . 

yi = XjXk\ r = 2a;ia;2a;3 = 2 1/2/1 2/21/3; Xi = . (83) 

D{xi,X2X3) yi 



IJ-i = 



Ai + A2 + A3 



-A,: 




exp 



^1 ^0 



'^1=1 



^ SiXjXk I = ;^ ( n / '^y^y'i^ ^ i-Siyi) I ( X] — 



=1 yi 



(84) 
(85) 
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Let us now use Mellin integral, summing up multinomial expansion 

_d 

(tl) ^ =^J D^f{n..)yT;Dn^^,.,= '^-.,^... (86) 

with integration contours Cr encircling poles of T{lij) at w = 0, —1, .... Integral is, in fact, symmetric with respect to the 
choice of two integration variables out of 101^102,103- 

/• 3 I* 00 f 3 

G{suS2,S3)=gi23f Dnl[r{io,) dy^y'^'+'^^-'e-'^y^ ^gu3f Dnl[sr'^'-'''r{Lo,)r {fi, + to,) , (87) 

This expression defines analytic continuation to complex values of Si. The remaining steps follow meromorphization of the 
2-point function, independently for each variable Si . We have for the corresponding Taylor coefficient in terms of : 

2m 

and finally, for the meromorphized triple vertex (with jii defined in ([M|) ): 

r(si,S2,S3) =M[G,Q,Sl,S2,S3]g(si)Q(s2)g(s3) = 7ni ,«„n3 4 ' ^3 % (89) 

ni,n2,n3>0 

/ duji f duj2 T-f , , , , d , , 

Inun^.ns = 9123 f f 77^-^\\1^,n,+^,i+i^iT{iO^)T{^i^+UJ^)]i03 = --UJl-U!2, (90) 

This is our goal: triple vertex defined as an entire function of each variable. The matrix indexes are to be inserted in obvious 
places, we omit them for brevity. The q coefficients inside the integral decrease as square of Gamma function of its arguments, 
providing decrease of the integral at large rii. This double contour integral can be reduced to double expansion over residues at 
LOi = —liU2 = —I2 of gamma functions T[loi)T{lo2) in left semiplanes 



7ni,n2,n3 = 5123 2_/ r(a;3) || gi,„^+p, T (^i + Wj) , = -^i , a;2 = -?2 , ^^3 = ^ + ^1 + ^2 ■ (91) 
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In particular, in conformal limit for q we obtain infinite sum of ratios of Gamma functions: 



oo 



3 



r {^i + uji) 



7, 



,0 

ni,n2.na 



5123 

hM=0 



■r(c.3)n 



i=l 



r(l + ni + /ij + uJi)T (1 + rii + 7i (A) + /ij + uJi) 



(92) 



In general case of running coupling constant, the proper way of meromorphization of the vertex function in perturbation 
expansion is to represent Feynman graphs as Mellin transforms |10j and meromorphize the powers of external variables Si as we 
did above, using g^. 



6 Spins and Operator Mixing 

So far we ignored the operator mixing, and studies idealized version of planar QFT, without spins. Let us now take both of these 
effects into consideration. As it was suggested in old papers [4] and reiterated recently [1] the proper framework, automatically 
preserving planar unitarity-analyticity is given by large order limit of matrix Fade approximant. We are going to refer to this 
limit as Matrix Meromorphization. Nothing changes in the general formulas of the Introduction, as long as we treat G, P, Q 
as infinite matrices in Hilbert space of composite fields made of quarks and gluons. It will become the Hilbert space of free 
composite particles of our planar QFT. 

Let us now go into details of this Matrix Meromorphization. The matrix G is acting in space of irreducible tensors of 
space-time symmetry group 0{d) (we are working in Euclidean space so far, but we do not assume conformal symmetry). There 
are tensor indexes for each of two operators averaged in G = (O1O2) • The tensor G is made of products of e^,yAp...j (^^i^ and 
kfj. = J^ with scalar functions of t — —k^ in front of these invariant tensors. In general, the ranks ni,n2 of Oi and O2 are 
different, so that G has rii + 712 indexes. For the same parity of Oi, O2 there will be no e^^Ap... tensors. 

As it was discussed at length in pj the Fade approximant. is nothing but a continued fraction summing up Taylor expansion 
near t = — in deep Euclidean region. In case there are tensor indexes, we can still treat G as analytic function of t with 
fixed unit vector k. There are some conspiracy relations for these scalar functions in front of jproducts of k which are needed 
to remove kinematical singularities at p = of the unit vector. In case of even number 2m of k factors (same parity of Oi, O2 
) the scalar function in front must vanish as otherwise it should have extra factor of ^/t. In order to avoid the kinematical 
singularities we must single out the factor i" and build Matrix Fade approximant. for t^"-G . It will have these kinematical 
poles at i = which will now cancel by t". 

The Fade equation for Q become matrix equation in this space of invariant tensors. 
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where L = M + 1,...,M + N and the contour R encloses the positive real axis clockwise (i.e. goes backwards from +00 to 
along the lower side of the cut, then forward from to +00 along the upper side of the cut. 

The Pade equations determine P(t), Q(t) up to arbitrary right multiplication by a tensor independent of t. This tensor can 
depend on unit vector k. This gauge invariance does not affect the matrix product P(t) • Q~^{t), which is our matrix Fade 
approximant. The problem of solution of Pade equations involve at some stage the problem of infinite matrix inversion for Q in 
this Hilbert space. 

In the leading conformal approximation, however, this problem of infinite matrix inversion dramatically simplifies, because 
the basis is known where the G matrix is block diagonal. This is the basis of irreducible conformal tensors. Corresponding 
2-point functions G{t) arc diagonal in all conformal quantum numbers, i.e. 0{d) quantum numbers plus scaling dimensions. 
In particular, in case of two symmetric traceless tensors of the same dimension A and the same rank n in coordinate space 
(different ranks or different dimensions do not correlate) 

/ n \ 

n-A 



G{x) = n (^'V-. - 2x^.a;.,) (x^) " ; (93) 



sym 



where sym denotes symmetrization and subtraction of traces for {/Xj} indexes as well as {vi} indexes. In momentum space 
(with 7 = A - d/2) 



G(fc) = a(n + 7) \J[{-dX,., +2d^,d^,) \ (fc2)"+^; (94) 

/ sym 

After all differentiations and symmetrizations we get 

n 

G(fc) = (fc2)^5^5;T;; (95) 
with complete set of invariant tensors T; and calculable coefficients gi in front (see Appendix2). Explicit form of these tensors 

(96) 

(97) 
(98) 
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Let us now multiply G{k) by Q{k). First, we expand Q{k) in the same set 

n 

Q(fc)= ^T„g(™)(t);t = -fc2 (99) 

m=0 

By construction these invariant tensors satisfy algebraic relation, with some Clebsch coefficients 

n 

TrT„ = 5]C/„T,- (100) 

3=0 

Note that Tr Tj • Tj ^ 6ij so that matrix inversion is needed to compute Clebsch coefficients from this relation. 

Ti = TvTi, 
T- ■ = TV T ■ • T • 
'-^ilm — Tr Ti • • Tm . 

In particular 

Tiio = Til, 

These matrices Tij,Tiim are independent of k in virtue of 0{d) invariance. We shall also use these Clebsch coefficients with lower 
indexes 

Cilm = TijCf^ = Tii^. (101) 
which are easier to compute. Making use of this algebra we get the set of equations (with 6 = 7 — n) 

n 

Aim = Y.9lTilm- (103) 
1=0 

We already know how to solve such equations by means the Greens function from Appendix 1. We write (with A = 1) 
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i-s)' 









Jr 27^^ 






i-tr 




jR 27rz 




Jr 27rz 



(104) 



From this equation (in the limit when M, N oo, = 2MN/A fixed, in units of R), we derive 

" m=0 r-=0 V ' V ' 

with some yet undetermined factor Xj . Now, expanding the entire functions in convergent series 

oo 

g(™)(s) = ^gM(-s)« (107) 



coefRcients qr"''^ as functions of r: 



and comparing coefBcients in front of {—ty and using F-^ir^u) = Sr,u we get finite set of n Unear algebraic equations for n 

n „ 
m=0 \ / V / 



, , (A-UX) 

9r"^ = T^7 TVF? — TT,m = 0,..n (109) 

r(r + l)r(r + 7-n+l) ^ ^ 

The constant vector X remains arbitrary here. This reflects the gauge invariance of right multiplication of Q(fc),P(A;) by 
arbitrary matrix W(A;). In our case this gauge matrix has a form 

n 

W(A;)=5^X^T,-. (110) 

The simplest choice would be 

Xi = Aio (111) 
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so that 



o(™) ^ (112) 

r(r+ l)r(r + 7-n+ 1)' ^ ' 



= 6,M-t)-'^Iy-n (2V^) (113) 
Let us now discuss the numerator of Matrix Meromorphized 2-point function: 

p(^)_ f ds(-s)-«G(s) . Q(s) 



Ct 27ri(s - t) 

Expanding in power series in k'^ — —t and integrating over s we get 



V{k) = T™P(")(-fc'), (114) 

{-ty 



^r(r + l)r(r-7 + n + l) V / 

The resulting meromorphized function reads 

M[G,Q,fc] = ,1^. ^ E-g^ n^-fc- n fc'^M.'^. (116) 

The mass spectrum m„ defined by roots of Bessel function 

J^_„_|(2mi?)=0. (117) 

For typical tensor family in QCD with some number of covariant derivatives A — n does not depend on n in zcroth order of 
perturbation theory. This dependence comes only from anomalous dimension as a function of running coupling. 

Above relations generalize the meromorphization equations for spin zero we have derived in the Introduction. The formulas 
(|17p do not rely upon conformal symmetry. As long as we use the Mellin transform for each scalar function in G{k) 

(fc2)-"G(fc) ^Y.^^f v^-FM htr , (118) 



25 



we get the following representation for P(fc), Q(fc) 



Q(fc)=To^g.fc2'- (119) 



r=0 



p{k) = Y,'rjJ2p'rk'^ (120) 

j=0 r=0 

4 = ti (121) 



m=0 



Equation for Q(fc) will change in case of running coupling constant. In the same way as we did it for spin zero case we can 
use the perturbation expansion in dimensional regularization in matrix form and relate terms of expansion to the corrections in 
general Pade equations for Q(fc). 
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8 Appendix 1. Missing Chapter of Pade Books 

Let us rederive here Greens function of Pade equations [J in modern notations and typeset. We express derivatives of our 
function Gq ~ (— t)''at t — —A as Cauchy integrals and we arrive at the following set of Pade equations (we set the normalization 
point A = 1 to simplify formulas) 



where L = M + 1,...,M + N and the contour R encloses the positive real axis clockwise (i.e. goes backwards from +00 to 
along the lower side of the cut, then forward from to +c» along the upper side of the cut. 

The discontinuity of {—ty along this cut equals to sin (tti/) which explains the factors in denominator. We are looking for 
the Greens function K{t, s) of Pade equations, which must satisfy inhomogeneous equations with proper right side: 



The implied extra condition is that K^{t, s) must be iV-th degree polynomial in t with s-dependent coefficients, so that we 
have linear system of integral equations for these coefficients as functions of s. In the same way, Qait) is iV-th degree polynomial 
in t with constant coefficients. 



th0906.4240 




(122) 




(123) 
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The general Pade equation, with arbitrary G(s) = {—sY{l + g{s)) reads 

ds-^^^{l + gis))Qis) = 0. (124) 

where Q{s) is A^-th degree polynomiaL Replacing by the left side of I123i this equation can be expressed as an 

integral equation using K, Qq: 

ds 

— K,{t, s) (1 + g{s)) Q{s) = Qo(i). (125) 
I R 27r« 

Let us check that the solution for K is given by the following Mellin-Barnes integral 

._{~iY'!if dz f dz' f{z){i + ty 1 

" fc 2^ fc 2^ f{z'){l + s)-'+i z - z' ' ^ ^ 

where 

/(.)= r(M + iv + i .)r(^z) 

■'^ r(Af + i-z^-z)r(iv + i-z)' ^ ^ 

and contour C encloses the poles of f{z) at z = 0, 1, ..N while contour C" encloses the zeroes of f{z') at z' = Af + 1 — v, M + 

2 — V, ... + GO. 

The same function /(z) determines the solution Qo of the homogeneous equation 

Qo(^)(xj^|i/(z)(l + ^)^ (128) 

Let us check this statement first. We obtain the following integral in (|122p 

—f(z) I = / —f(z)B (L-iy-z,iy+l) (129) 

c2^z^^ Vfl2isin(7ri.)(i + i)^+i-- Jc ' V , ; V 7 

The function /(z) was chosen in such a way that the integrand here reduces to the rational function 

M+N L-1 N 

f{z)B{L-v- z,v+l)o^ W{k-z) Jl {l-v - z)W^{n- z)-^. (130) 

k=L+l l=M+l n=0 

This function does not have poles outside integration contour C and decreases as z~^ at infinity, therefore integral is equal to 
zero. On the other hand, /(z) by itself has poles at z = 0, 1, ...N which makes Qq a polynomial of N — th degree. So, this is a 
solution of Pade equations. 
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When we substitute the Anzatz (|126p back into (|123p we first integrate over t and we get this time 

fiz)B{L-iy-z,iy+l) 



z — z 



(131) 



which has exactly one pole at z — z' outside integration contour C and decreases as z~^ at infinity. Therefore, the z 
integration reduces to the residue at this pole. The factors f{z) and f{z') cancel among themselves and we get the standard 
integral, calculable by taking residues at z' = L — v — n. Summing up resulting binomial expansion we finally get 

, / dz' B{L-V~Z',V + 1) (-l)-"r(^ + 1) ,.-n-L-l_ (1 + ^)'^ M X ^ i-'^y (,on^ 

^ ' Jc,27Ti (l + s)-'+i n!r(i.-n + l) ' "(1 + 5)^+1^ + 5 7 "(l + s)^+i' ^ ' 

which is the RHS of ([m)) . 

Now, let both Af , go to oo. In this limit (assuming z ^ MN) 

f{z) ^ i-zr-' exp (-^) , zo - 2MN. (133) 

By rescaling variables (remember, there were factors of A for each t, s variable, so now we switch to the units where R ~ I): 

I I t s 

Z^ZoZ,Z ZqZ ,t ^ — , s — 

zo zo 

we arrive at the integrals 

dz r (iz'exp(i-i+tz-sz')(-z) 
sin(7rz^) Jfj 2TTi Jqi 2iTi {z — z') (z') 

Qo(0 / ^{-zf-' exp (-- +tz)=Y, J~T^,^ (135) 
Jc 27ri \ z J n\T{n + v + I) 

In this paper we use the expansion coefficient Ki,^m{s) of Ki,{t,s) in front of (— i)™, integrated over s with weight (— s)", 
corresponding to Mellin transform 

1,2^ "'"'^ '^ ' " sm{7ru)m\T i-u;) Jc 2Tn fc' [z - z') {zr+^ ^ ' 

This integral simplifies by the following change of variables 

z^- z'^ — 
p' q 
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and adding one more integration 

Hdxcf ^(f ™(-q)-+"-iexp(-(l + x)(p + g)), 

Jo Jr^ 27ri Jj^ 2m 

where contour R encircle positive axis clockwise, as before. These integrals reduce to Gamma functions: 



^ ^(--Z) exp + .)) = 

^(-p)— exp(-p(l + x)) = + 

^ Am i (to + j/j 

' 1 







w — m 



(137) 



Collecting all factors we arrive at expression (P7| 

£ |^^-.™(«)(-^)" = /^-K^)' (138) 

Taking the limit at integer w — n > 0, 

/ :^if.,™(s)(-.s)" = (5„™. (139) 
jR 27ri 

This normalization provides that at 5(5) = our integral equation (jl25p has solution (3(s) = (9o(s). 



9 Appendix 2. Conformal 2-point function in momentum space 



The momentum space 2-point function of two symmetric traceless tensors of rank n has the form (|94p . 

In order to compute these gm one may use generating function by multiplying this equation by products of complex light-like 
vectors ^, r\. 

All the trace terms disappear in virtue of = 77^ = and we are left with the following polynomial 

n 

a (n + 7) (-(c^)9' + 2 m mT i^y^" = (k'r E 3- (^^)™ {^>'T (^^)"""' (i^o) 

m— 
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Let us introduce the polynomial F; {x) as follows (we keep constant parameter 6 = 7 + n) 

Fa.) = aib){k^y-''[-d^+'-^^)\kr, (141) 
i^k) (nk) , , 

n 

F„(a;) = QmX^. (143) 

m— 

Shifting / by 1 we find recurrent equation 

F,,. (.) = [er-^ {-d^ + ^i^) (.^)-' (.) , (144) 
Fo{x) = a{b). (145) 
We use identities (with d being gradient with respect to fc ) 

(r?a)x=M(i_2x), 
(C5)x=M(i_2ar), 
((r?a)x)((C5)x) = '=(^~^'^^' 



(dx) 



{jri) ((1 - 2x) - 2x{\ - 2x) - 2x(l - 2x)) _ {$,r}){l - 2x){l - Ax) 
2 + + 4x2 + 2a; - 4a;2 - 4x2 2x(l - 2x) 



fc2 fc2 

(fca) X = 0, 
2 

Q^x = — (1 — dx) . 
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We find 

(e)'-' Fi {x) = d (2(6 - l)k (fc2)'-'-i Fi (x) + (fc^)"-' dxFl (x)) 

= 2{b - l){d + 2(6 - ; - 1)) Fi (x) + 0+ 

+ (e)"-' d'^xFf (x) + (fe2)"-' {;dxf F[' (x) 

= {k'^f"^'^ (2(& -l){d + 2{h-l- l))Fi {x) + 2(1- dx) Fl (x) + 2x{l - 2x)F/' {x)) (146) 



and 



' Fi{x) = m (2(6- 



m {k^t 

I (2(6- 



F,(x) + (fc2)'^'((r?a)x)^^/(x)) 

-0+4(6-0(&-^-l)a;)Fi(x) \ 
+2(6-0a;(l-2x)F/(x) 
+2(6- Z)x(l - 2x)F[ \x) 
+(l-2a;)(l-4x)i^/(x) 

+a;(l - 2a;)2F/' (x) / 



Combining these terms we get 



Fi+i (x) = (-2 (6 - (26 + d - 2Z - 4) - 8 (6 - (^ - & + 1) x) Fi (x) (147) 

+ (2 (46 + d - 4? - 6) X - 16 (6 - ^ - 1) x^) F/ (x) (148) 

+ (2x - 1)4x2f/' (x) (149) 

Finally, the recurrent equation reads 
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(X) 


- ( 40 ^ 

- + 


AiX + + Hi X) X — 


^? 


= -2{b- 


-l){2b + d-2l-4); 


Al 


= -8(6- 


~l){l-b+l), 


Bf 


= 2 (46 + 


d-4l-6); 


Bl 


= -16(6 


-l-l), 


CO 


= -4; 






= 8; 




Fo{x) 


= a{b). 





C"" + C'x)x^^j Fi{x); (150) 

(151) 
(152) 
(153) 
(154) 
(155) 
(156) 
(157) 

Above differential equations provide recurrent relations which allow one to find expansion coefficients fm one after another, 
for I = 0,1, ...n. 

= M„,„_i(/)£Li + Mrn,m{l)f^\m = 0, ...I + 1, (158) 
/I'l - //|\ = 0; (159) 
Mm,mil) = ^? + mB^ + m{m - 1)C°, (160) 
Mm,m-i{l) = Aj + (m - l)Bl + (m - l)(m - 2)C^ (161) 

This is multiplication by bi-diagonal Z x (/ — 1) matrix M(Z) with diagonal elements Mm,m{0 and sub-diagonal elements 
Mm,m-i{l) , other elements being equal to zero. So, the solution of this recurrent equation can be written as a matrix product 
for vector g with components Qm 

g = f(") = M(n - l)M(n - 2)...M(0)f(°\ (162) 

= a{b)SmO, (163) 
6 = 7 + n. (164) 
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